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The quarkonium hadroproduction in association with a photon at high energies provides a probe 
of the dynamics of the strong interactions as it is dependent on the nuclear gluon distribution. 
Therefore, it could be used to constrain the behavior of the nuclear gluon distribution in proton- 
nucleus and nucleus-nucleus collisions. Such processes are useful to single out the magnitude of 
the shadowing/antishadowing effects in the nuclear parton densities. In this work we investigate 
the influence of nuclear effects in the production of J/tjj +7 and T + 7 and estimate the transverse 
momentum dependence of the nuclear modification factors. The theoretical framework considered in 
the J/tp (T) production associated with a direct photon at the hadron collider is the non-relativistic 
QCD (NRQCD) factorization formalism. 

PACS numbers: 12.38.Bx, 13.25Gv,13.60.Le, 13.85.Qk; 12.38.-t 
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I. INTRODUCTION 

In recent years, quarkonium hadroproduction has be- 
come the subject of intense theoretical and experimental 
investigation. The main reason is that production and 
decays of heavy quarkonia have been an ideal labora- 
tory to investigate Quantum Chromodynamics (QCD). 
Their large masses provide a hard scale which allows us 
to use perturbative QCD techniques. Basically, there 
exist three distinct formalisms for quarkonium produc- 
tion. The simplest one is the color evaporation model 
(CEM) where the hadronization of the QQ pairs into 
quarkonia is assumed to be dominated by long-distance 
fluctuations of gluon fields which motivates a statistical 
treatment of color. In the color-singlet model (CSM) |2j, 
quarkonium is viewed as a non-relativistic color-singlet 
bound state of a QQ pair with definite angular momen- 
tum quantum numbers. Production rates for a quarko- 
nium state are computed by calculating the production 
of a heavy quark pair which is constrained to be in a 
color-singlet state and have the same angular momentum 
quantum numbers as the physical quarkonium. In the 
non-relativistic QCD (NRQCD) factorization formalism 
3], non-perturbative aspects of quarkonium production 
are organized in an expansion in powers of v, the rela- 
tive velocity of the QQ in quarkonia. For production of 
S-wave quarkonia, the results of the CSM are recovered 
in the limit of v — > 0. In addition, in the NRQCD for- 
malism new quarkonium production mechanisms are now 
possible since it is no longer required that the QQ pro- 
duced in the short-distance process have the same color 
and angular momentum quantum numbers as the quarko- 
nium state. Recently, several computations of the QCD 
corrections to the inclusive quarkonium hadroproduction 
processes have shed some light on the robustness and/or 
deficiencies of those models. Comprehensive reviews of 
recent developments in the theory of quarkonium pro- 
duction can be found in Refs. 0, [f|. 

Beside the studies of inclusive production, efforts are 



being made to obtain improved theoretical predictions for 
complementary observables to the inclusive yield, such 
the hadroproduction of J / ip and T in association with a 
photon. In the framework of CSM, the associated pro- 
duction of J ftp + 7 at a hadron collider was first proposed 
as a good channel to investigate the gluon distribution in 
the proton with a relatively clean signal @. In Ref. 0], 
such a process at the Tevatron energy has been consid- 
ered in the CSM at LO, and the results show that the 
contribution from the gluon fusion sub-process is domi- 
nant over that from the fragmentation process (the same 
occurs at the LHC Q). The color-octet contributions 
were investigated in Ref. @ and it was found they are 
dominant in the large px region. Recently, in Ref. [Toj 
the effect of the NLO QCD corrections to J ftp (T) + 7 
hadroproduction at the LHC has been investigated. In 
Ref. the real next-to-next-to-leading (NNLO) order 
QCD contribution to hadroproduction of a J/ip (T) + 7 
via color singlet transitions for the inclusive case has been 
addressed. In this Letter, we examine the production of 
associated J/ip + 7 and T + 7 at large pt within the 
NRQCD approach in proton-nucleus and nucleus-nucleus 
collisions. Such processes are relatively clean because the 
produced large pt quarkonium is easy to detect through 
its leptonic decay modes and the quarkonium's large pt 
is balanced by the associated high energy photon. At 
the LHC energy, the leading contributions are dominated 
by gluon induced hard processes and then the quarko- 
nium production associated with a direct photon will be 
strongly dependent on the nuclear gluon distribution. 

Our goal is to use the J/ip + 7 and T + 7 processes as 
auxiliary observables to constrain the nuclear gluon dis- 
tribution. This is motivated by similar investigations on 
inclusive heavy quark, quarkonium and prompt photon 
production in central proton-nucleus and nucleus-nucleus 
collisions (See e.g. Refs. [1343 )■ One of the nuclear ef- 
fects which is expected to modify the behavior of gluon 
distribution is the nuclear shadowing. This effect has 
been observed in the nuclear structure functions by dif- 
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ferent experimental collaborations [13, [13] in the study 
of the deep inelastic lepton scattering (DIS) off nuclei. 
The modifications on F%(x, Q 2 ) depend on the parton 
momentum fraction x. While for momentum fractions 
x < O.f (shadowing region) and 0.3 < x < 0.7 (EMC 
region), a depletion is observed in the nuclear structure 
functions, in the intermediate region (O.f < x < 0.3) it is 
verified an enhancement known as antishadowing. These 
experimental results strongly constrain the behavior of 
the nuclear quark distributions, whereas the the nuclear 
gluon distribution is still an open question due to the 
scarce experimental data in the small- a; region and/or 
for observables strongly dependent on the nuclear gluon 
distribution. 

This Letter is organized as follows. In next section, 
we summarize the main formulas concerning the process 
J/ip + 7 in the NRQCD formalism and define the nu- 
clear modification factors for proton-nucleus and nucleus- 
nucleus collisions, R p a and Raa, respectively. In last 
section we present the numerical results considering the 
more recent nuclear parton parameterizations and esti- 
mating the transverse momentum dependence of the nu- 
clear modification factors at the LHC energies. 

II. QUARKONIUM PRODUCTION 
ASSOCIATED WITH A DIRECT PHOTON 

As long as J/ip + 7 is produced at small longitudi- 
nal momentum fraction, xp <C I , the gluon fusion chan- 
nel dominates over the qq annihilation process. There- 
fore, at high energies and at leading order (LO), the pro- 
cess g + g — > J/ip + 7 contributes at the partonic level 
with six Feynman diagrams, which is similar to that of 
g + g — > J/ip + g in inclusive J j ip hadroproduction. The 
signal we focus on is the production of a J/ip and an 
isolated photon produced back-to-back, with their trans- 
verse momenta balanced. The LO cross section is ob- 
tained by convoluting the partonic cross section with the 
parton distribution function (PDF), g(x,fip), in the pro- 
ton, where fip is the factorization scale. At NLO expan- 
sion on a s , there are one virtual correction and three 
real corrections processes, as shown in Ref. [13]. In 
the NRQCD formalism, the contributing subprocesses 
are q + q -> 2S+1 Lj + 7 and g + g -t 2S+1 Lj + 7. 
The Fock-components that contribute to J/ip produc- 
tion are the color-sing let 3 S[ 1] state and the color-octet 
states 3 s[ 8] , l sf ] and 3 p£I j2 . The color-singlet 3 5i state 
contributes at 0(1) but the color-octet channels all con- 
tribute higher orders in v (the relative velocity between 
the heavy quarks). 

We are interested here in the quarkonium production 
in a nuclear medium. In order to get the J /ip + 7 yield in 
pA and AA collisions, a shadowing-correction factor has 
to be applied to the J/ip yield obtained from the simple 
superposition of the equivalent number of pp collisions. 
This shadowing factor can be expressed in terms of the 
ratios Rf of the nuclear parton distribution functions in 



a nucleon of a nucleus A to the PDF in the free nucleon. 
Most of shadowing models provide the nuclear ratios at a 
given value of Qq and then evoluted through the DGLAP 
evolution equations [2l[ to LO accuracy. Only very re- 
cently, the nuclear PDFs have been available at NLO ac- 
curacy. Therefore, in what follows we will consider a LO 
calculation for the nuclear modification factors in order 
to be consistent with the limitation of shadowing models. 
In this respect our analysis cannot be considered as fully 
NLO and should be updated once the NLO calculation in 
NRQCD approach is available for the nucleon case. On 
the other hand, this is not an important limitation to our 
results as we will discuss later on. 

In order to obtain the transverse momentum (pr) dis- 
tribution for the process .<? + .<?—> J/ip + j, we express the 
differential cross section as 

(Per f 2 2 4:X!X 2 Pt da 

— — = / dx 1 g A {x u n F )g B {x2,H F )- — --tj, (1) 

dydpr J 2xi - x T eV dt 

where we have defined xt = Imr/y/s, with y/s be- 
ing the center of mass energy of the AB system and 

rriT = \Jpt + rri^ being the transverse mass of outgoing 

J/ip. The gluon distribution, gA/B(x, Q 2 ), in the hadron 
A/B is evaluated at factorization scale fiF- The common 
transverse momentum of the outgoing particles is px and 
y is the rapidity of outgoing J j ip having mass . The 
variables 27 and x 2 are the momentum fractions of the 
partons, where M 2 / s < x\ < 1 (M is the invariant mass 
of J/ip + 7 system) and x 2 can be written in terms of 
other variables as 

Xl x T e-y - 2t m% 
x 2 = -z = — — , with r = — v -. 2) 

2x\ — xre y s 

In the NRQCD formalism, the cross section for the 
production of a quarkonium state H is written as <j(H) = 
J2 n c n(0\On\0), where the short-distance coefficients c n 
are computable in perturbation theory. The (0|O^|0) are 
matrix elements of NRQCD operators of the form 

(0\O% |0) =J2J2^n\H(\)+X)(H(\)+X\K n \0). (3) 

X A 

The bilinear in heavy quark fields which cre- 

ates a QQ pair in a state with definite color and an- 
gular momentum quantum numbers. Hereafter, we will 
use a shorthand notation in which the matrix elements 
are given as (0% s s ) ( 2S+1 Lj)) . The angular momentum 

quantum numbers of the QQ produced in the short- 
distance process are given in standard spectroscopic no- 
tation, and the subscript refers to the color configuration 
of the QQ: 1 for a color singlet and 8 for a color octet. 

The parton level differential cross sections relevant for 
hadroproduction of J/ip + 7, including both color-singlet 
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FIG. 1: Transverse momentum dependence of the nuclear modification factor R p a in J/ip + "f production in central and forward 
rapidities at the LHC (y/a = 8.8 TeV), for distinct nPDFs. 



and color-octet contributions are given below 0, [22| : 
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InEqs. (gj) we have defined <Jo = ir 2 e 2 aa 2 m c js 2 (with 
charm quark mass m c = 1.5 GeV) and si = s — Am 2 , 
t\ = t — Am 2 , and u± = u — Am 2 . In these formulae, 
s, i, and u are the Mandelstam variables, which can be 
written as 



X1X2S, t = m 



x\\/ smT& 



In our numerical calculations the one loop expression 
for the running coupling, a s (^R), with Aqcd — 0.2 GeV 
and n/ = 4 is considered. The (renormalization and fac- 
torization) scale for the strong coupling and for the eval- 
uation of PDFs is fj,p = fi^ — (pj, + m\), where m^ is 
the J/ip mass. For numerical values of the NRQCD ma- 
trix elements we have used those from Ref. [23], which 

are (units of GeV 3 ): (0( /lp ( 3 Si)) = 1.16, (O^fSi)) = 

1.19 x 10- 2 , (Og^^So)) = (Og^V^o)}/™ 2 - 0.01. Wc 
have checked that using another set of color octet matrix 
elements, taken from Q, our results do not change con- 
siderably, since we only calculate the ratio between the 
production cross sections, and the dependency on the 
color octet matrix elements cancels out in those ratios. 
In what follows we estimate the differential cross sections 
for central (y = 0) and forward (y = 2) rapidities for the 
pp, pA and AA collisions in order to compute the nu- 
clear modification factors. From Eq. (p}, it implies that 



the differential cross section at small values of the J/ ip 
transverse momentum in pPb (and PbPb) collisions at 
the LHC is determined by the behavior of the nuclear 
gluon distribution at X2 ^ 10~ 4 . It is important to em- 
phasize that smaller values of x contribute at J/iJ) + 7 
production in the forward rapidity region which can be 
measured, e.g., with the CMS and ALICE experiments 
at LHC. 

The main input in the calculations of J/ , 0+7 cross sec- 
tions is the nuclear parton distribution function (nPDF). 
Over last years several groups has proposed parameter- 
izations for the nPDFs 24-28], which are based on dif- 
ferent assumptions and techniques to perform a global 
fit of different sets of data using the DGLAP evolution 
equations. These parameterizations predict very distinct 
magnitudes for the nuclear effects. For larger values of 
x, the EKS and the EPS nPDFs show antishadowing, 
while this effect is absent for the HKN and DS parame- 
terizations in the x < 10 _1 domain. While the nuclear 
shadowing is moderate for DS and HKN parameteriza- 
tions and somewhat bigger for EKS one, the EPS predic- 
tion has a much stronger suppression compared with the 
other parameterizations. For smaller x around x ~ 10~ 5 , 
while DS and HKN parameterizations have about 20% 
suppression and EKS one have about 40% suppression, 
for the EPS parameterization this effect goes to almost 
80% suppression in the nuclear gluon compared with the 
A scaled gluon content in the proton. 

In what follows we calculate the quarkonium produc- 
tion in association with a direct photon in pA and AA 
collisions, considering the nuclear parton distributions 
discussed above. We then estimate the nuclear modi- 
fication factor for these predictions. These factors could 
be measured at the LHC and we analyze the transverse 
momentum dependence of them. Since the full NLO cal- 
culation of J/fy(T) +7 including the COM contributions 
is not available yet, and for the sake of simplicity, the cal- 
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FIG. 2: Transverse momentum dependence of the nuclear modification factor Raa in J/tp + y production in central and forward 
rapidities at the LHC (y/s = 5.5 TeV), for distinct nPDFs. 



culation of the cross sections is here done at LO accuracy. 
In order of to be consistent, we consider only the LO ver- 
sion of the nPDFs employed (besides, the EKS nPDF is 
only evolved to leading order) . As we are calculating ra- 
tios between cross sections, common uncertainties on the 
normalization of the pA, AA and pp cross sections, e.g. 
due to higher order contributions, are expected to cancel 
out in the nuclear modification factors. 

In Fig. Q] we present our estimates for the transverse 
momentum dependence of the ratio R p a for central and 
forward rapidities, defined by 



R 



pA 



_ da(pA) 



dyd 2 pj 



da(pp) 



(5) 



in pA collisions at the LHC (y/s = 8.8 TeV). The results 
show distinct behaviors of RpA for different nPDF's. The 
difference between central and forward rapidities for R p a 
comes from the kinematical xi range probed in the two 
cases. While for central rapidities the values of x^ are 
ever larger than 10 -3 , in the forward case the minimum 
value could be 10~ 4 , increasing with pr- Consequently, 
the effects in the nuclear gluon distribution which con- 
tribute for the quarkonium production are different in the 
two rapidities. In the pt interval shown in Fig. [TJ anti- 
shadowing effect is not observed as X2 < 0.1. We have 
that the nuclear factor is substantially suppressed in the 
EPS (EKS) case, going down an 0.8 - 0.7 at p T ~ 10 
GeV, while in the DS (HKN) case it is almost flat and 
equal to 0.95 (0.85) in the full pr range. Moreover, dif- 
ferently from the DS and HKN predictions, the EKS and 
EPS parameterizations lead to a strong transverse mo- 
mentum dependence. Consequently, the determination 
of the magnitude and px dependence of the this nuclear 
modification factor at the LHC could be useful to deter- 
mine the properties of the shadowing in the gluon distri- 
bution. 

The production of J/ip + 7 can also be studied in the 



collision of two heavy nuclei. In what follows we present 
our estimates for the transverse momentum dependence 
of the nuclear modification factor Raa for central and 
forward rapidities denned by 



Raa 



da {A A) 



dyd 2 p 7 



/A'- 



da(pp) 



dyd 2 p 7 



(6) 



in A A collisions at the LHC (y/s = 5.5 TeV). In this case 
the nuclear effects are amplified by the presence of two 
nuclei in the initial state of the collision. Similarly to 
the pA case, we can see in Fig. [2] that the factor Raa is 
strongly modified by the shadowing effects. For low px 
the suppression is stronger than in the pA case. We can 
see an anti-shadowing effect appearing in central rapidity 
for EPS and EKS nPDFs in the region of larger pr — 
30 GeV. Our results indicate that J/ip + 7 production 
with pt < 20 GeV can be used to determine the gluon 
shadowing effect. 

The production of T(15*) +7 can be obtained from the 
expression Q above, by replacing the charm mass and 
charge by the bottom ones, mb = 4.7 GeV, eb, the J/ip 
mass by the T(15) mass, and by using the correspond- 
ing color octet matrix elements. We use the values taken 
from H|, namely (units of GeV 3 ): (Oj ( 3 S*i)) = 10.9, 
(O s T ( 3 Si)> = 0.02, (Oj^So)) = 0.136 and (OJCPq)) = 
0. We notice that using the alternative set for the 
color octet matrix elements [Hj]: (Oj( 3 5i)) = 0.03, 
(Oj^So)) = and 5(O|'( 1 P ))/m2 = 0.139, our results 
are not significantly changed (there is only a very tiny 
reduction in R for larger pr)- 

In Fig. [3] we present our estimates for the transverse 
momentum dependence of the ratio R p a for central and 
forward rapidities, in pA collisions at the LHC (y/s = 8.8 
TeV) producing T(IS) + 7. In this case the pr range is 
extended to 100 GeV, since the distributions tend to be 
shifted to larger values of px due to the larger T mass, 
which makes larger values of x to be accessed, and sug- 
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FIG. 3: Transverse momentum dependence of the nuclear modification factor R p a in T(15) + 7 production in central and 
forward rapidities at the LHC (y/s = 8.8 TeV), for distinct nPDFs. 



gests the existence of antishadowing. As in the previous 
cases the results show distinct behaviors of R p a for dif- 
ferent nPDF's. For central rapidities, the EKS and EPS 
show antishadowing for larger px, while this effect is ab- 
sent for forward rapidities and for the other distributions. 
For lower px, the suppression is different for the different 
nPDFs, being stronger in the forward case. This effect 
could then be used to discriminate among the nuclear 
PDF's. 



The results for the production of T(1S) + 7 in A A 
collisions are shown in Fig. 2J where we present our re- 
sults for the transverse momentum dependence of the 
ratio Raa for central and forward rapidities at the LHC 
(\/s = 5.5 TeV). Here the nuclear effects are more pro- 
nounced than in the pA case. The EPS parameterization 
shows the steepest behavior, presenting shadowing (sup- 
pression) for low pt and antishadowing (enhancement) 
for larger px and central rapidities. The EKS has a sim- 
ilar although less pronounced behavior, and the HKN, 
which presented an almost flat behavior in the J/ip case, 
it now grows steeper. Thus, the HKN versus DS could 
be discriminated at low px and central rapidities, and at 
high px and forward rapidities. In the forward case, anti- 
shadowing behavior at high px tend to be less important 
for the EKS and EPS nPDF's as px grows. We conclude 
that shadowing and antishadowing effects could be tested 
in different px and rapidity regions and therefore help in 
discriminating among the different nuclear PDF's. 



Finally, we discuss qualitatively the results and com- 
ment on the limitations of present calculations. Let us 
first concentrate on the pA case where cold matter effects 
also play an essential role. At leading order, the differ- 
ential cross section is simply proportional to the product 



of gluon densities 

, d<J , (pA -> J/-0 + 7) oc gp(xi,Q 2 )g A (x 2 ,Q 2 ) 

CLX 1(2X2 

XS(XIX 2 S - 771^), 

where x% t 2 are the projectile and target-parton momen- 
tum fractions. In this kinematics, the nuclear modifi- 
cation factor reduces to R p a — Rg(x2)- In our analy- 
sis above we did not take into account the nuclear ab- 
sorption. In the framework of the probabilistic Glauber 
model, this effect refers to the probability for the pre- 
resonant QQ pair to survive to the propagation through 
the nuclear medium. Therefore, the J/ip may be sensi- 
tive to inelastic rescattering processes in a large nuclei, 
which spoil the simple relationship between R p a and R^. 
Assuming the factorization between the quarkonium pro- 
duction process and the subsequent possible J/ip inelastic 
interaction with nuclear matter, we can write the produc- 
tion cross section as 

^(p a -> J /^P + l) K S ahs (A, cr//,/,^) 

K^{pA^J^ + 1 ), 
dx 2 

where S^s denotes the probability for no interaction, or 
survival probability, of the meson with the nucleus target. 
In a Glauber model it depends on the J/ip — N inelastic 
cross section, (Jj/^nj and reads as (for large nucleus): 

S a bs * -r-^ / d 2 b [1 - exp (-T A (b) aj^ N )] ,(7) 

A<J J/ipN J 

with the thickness function T A (b). As ctj/^n is n °t well 
constrained from data, this is additional uncertainty en- 
tering on the nuclear modification factor. In a rough 
estimation we have R p a — S a Y >s (A,aj/^ N )R^. The cor- 
responding expression for AA collision can be obtained 
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using similar methods. We quote Ref. [30( as an exam- 
ple where the magnitude of these corrections has been 
studied for the J/tp production in proton-nucleus and 
nucleus-nucleus collisions. 

Concerning nucleus-nucleus collisions, the situation is 
more complicated as final state effects can not be disre- 
garded. For instance, the processes of dissociation and 
recombination of cc pairs in the dense medium can be 
computed through the co-movers interaction model [3l[ , 
which also incorporates also the recombination mecha- 
nism }32j [. As mentioned for pA case, nuclear effects in 
nucleus-nucleus collisions are usually expressed through 
the so-called nuclear modification factor, R'^g(b), de- 
fined as the ratio of the J/ ip yield in AA and pp scaled by 
the number of binary nucleon-nucleon collisions, N co u(b). 
A similar factor can be defined in our case of J/\& + 7 
production, having in mind that the prompt photon is 
insensitive to nuclear matter effects. We have then for 
symmetric nuclei, 



N co n(b) dN^/dy 
J d 2 s (JAAjb) n(b, s) P sup (b, s) 
J d 2 sa AA {b)n(b,s) 



(8) 



where (JAA(b) = 1 — exp[— a pp A 2 TAA{b)] and T A A{b) = 
/ d 2 sTA(s)TA(b — s) is the nuclear overlapping function. 
One has that T A (b) is obtained from Woods-Saxon nu- 
clear densities, and 



,(M) = s s j%(b, s )s abs (b, s )s co (b, s ), 



(9) 



n(b,s) = <T pp A 2 T A {s)T A {b- s)/a AA {b) , (10) 

where the number of binary nucleon-nucleon collisions at 
impact parameter b is given by, N co n(b) = f d 2 sn(s,b). 
The three factors appearing in Eq. ©, S sh , S abs and 
S co , denote the effects of shadowing, nuclear absorption, 
and interaction with the co-moving matter, respectively. 

As referred above, the nuclear absorption is usually in- 
terpreted as suppression of J/ tp yield because of multiple 
scattering of a cc pair within the nuclear medium. At low 
energies the primordial spectrum of particles created in 
scattering off a nucleus is mainly altered by interactions 
with the nuclear matter they traverse on the way out to 
the detector and energy-momentum conservation (it has 
been shown that it becomes a minor effect at AA colli- 
sions at LHC). For nucleus- nucleus collisions these effects 
can be combined into the generalized suppression factor 



gabs 



[l-exp(-P(xi)acgAT A (b))] 



(3(x 1 )(3(x 2 )a 2 c5 A 2 T A (s)T A (b 
x [1 -aip(-l3(x 2 )a c£ AT A (b- 



-s) 

«))] 



(11) 



where Xi t 2 = (v^l — 4M 2 /s ± xp)/2, and (3(xi^) = 
(1 — e) + exj 2 determines both absorption and energy- 
momentum conservation. The parameters 7, e and tr cS 
can be adjusted to describe collider data. Secondly, co- 
herence effects will lead to nuclear shadowing for both 



soft and hard processes and therefore for the production 
of heavy flavor. Shadowing factor S s J]^{b, s) can be cal- 
culated within the Glauber-Gribov theory, and here it is 
already included in the nuclear PDF ratios. 

The final state effects can be addressed taking, for ex- 
ample, the co- movers interaction model 32[. Assuming a 
pure longitudinal expansion and boost invariance of the 
system, the rate equation which includes both dissocia- 
tion and recombination effects for the density of charmo- 
nium at a given production point at impact parameter s 
reads 



dNj/^(b,s,y) 



dr 



Oco 
T 



N c °{b 1 s 1 y)Nj / ^{b 1 s 1 y) 
N c (b,8,y)N g (b,a,yj\ , (12) 



where N co , Nj/^ and N c tc) is the density of comovers, 
J/ip and open charm, respectively, and a co is the inter- 
action cross section for both dissociation of charmonium 
with co-movers and regeneration of J/if) from cc pairs 
in the system averaged over the momentum distribution 
of the participants. It is the constant of proportionality 
for both the dissociation and recombination terms due 
to detailed balance Nj^N co = N C N 5 . The solution of 
Eq. (fl"2|) can be approximated by 

N co 



S co = cxp \-(T co [N c °(b,s,y)-N ah (b,s,y)] In 



N pp (0)_ 



N 



sh 



(dcT pP /dy) 2 
o^da^jdy 



N bin (b,s)S shad (b,s,y) 



(13) 



Details of the model can be found in [32| and the quan- 
tities in Eq. (TT5)) are all related to pp collisions at the 
corresponding energy and should be taken from exper- 
iment (the extrapolation for LHC is quite uncertain). 
The magnitude of recombination effect is controlled by 
the total charm cross section in pp collisions and then 
dissociation-recombination effects will be of crucial im- 
portance in PbPb collisions at LHC. 

At this stage is difficult to single out the size of each 
contribution at the LHC for the J/\& (T) + 7. For such 
photon-associated high-p^ quarkonium production, the 
final production cross section is modified by nuclear 
absorption, or final-state interaction between the pre- 
resonance state and the nuclear medium as discussed 
above. Notice that the Glauber absorption model is only 
valid for total production cross section which is domi- 
nated by low-pT quarkonium production. Here, we fo- 
cus on the transverse momentum dependence of the nu- 
clear modification factor due to gluon shadowing and it 
is timely to discuss the px dependence of the final-state 
interaction. To address this issue, we can do an educated 
guess by verifying what we learn from RHIC. The mod- 
ification factor R AA has been determined as a function 
of pt in the central rapidity region, analyzed both by the 
PHENIX [33] and STAR [H collaborations. The fac- 
tor increases with px and STAR data indicate a nuclear 
modification factor larger than one for pt ~ 8 GeV [34[ • 
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FIG. 4: Transverse momentum dependence of the nuclear modification factor Raa in T(15) + 7 production in central and 
forward rapidities at the LHC (y/s = 5.5 TeV), for distinct nPDFs. 



If this non-suppression of J/tp at large px is confirmed it 
does not seem to behave as the other hadrons, which are 
significantly suppressed in central collisions and for in- 
creasing pt- This behavior is closer to the one observed 
in prompt photons production. In the case of quarko- 
nium plus a prompt photon, we then would expect the 
same non-suppresion pattern, since the hard photon is 
not affected by nuclear matter effects. 

As a final comment, we notice that the present calcula- 
tion can be also extended for the RHIC kinematic region, 
where the qq annihilation processes could also contribute. 
For instance, a study has been done in Ref. [IH] for the 
inclusive J / ip production in pp collisions at RHIC within 
the PHENIX detector acceptance range using the color 
singlet and the color octet mechanisms based on NRQCD 
formalism. The calculations reproduce the RHIC data 
with respect to the px distribution, the rapidity distri- 
bution and the total cross section at y/s — 200 GeV. In 
Ref. [30J , the authors have considered the CSM approach 
(taking into account the s-channel cut contributions (36[) 
and cold nuclear matter effects to compute the J/ip nu- 
clear modification factor for dAu and AuAu/CuCu col- 
lisions at RHIC energies. The data description is quite 
satisfactory. Therefore, this indicates the robustness of 
the framework in describing the high energy data. Of 
course, we are aware of the limitations of present calcula- 
tion. The picture of a J/tp exactly recoiling back to back 
with the photon is simplified, since possible contributions 
of the color-octet channel would lead to a fragmentation 
of the gluon into the meson, which would not be exactly 
recoiling against the photon. The experimental detection 
could be an additional concern. For example, ALICE 37] 
can tag the muons in the forward region, where however 
there is no photon reconstruction, while it can detect the 



photons in the central region, where however the muon 
detection is not as good. In the case of CMS (38|, the 
thresholds in px both for muons and for the photon ap- 
pear quite low to ensure an efficient triggering. 

In summary, in this letter we have investigated the 
quarkonium production in association with a isolated 
photon in pA and A A collisions at the LHC, considering 
the NRQCD factorization formalism and some of the pa- 
rameterizations for the nuclear parton distributions avail- 
able in the literature. Our results show that the nuclear 
modification factors R p a and Raa are useful observables 
to discriminate among the different parton distributions 
in the nuclear medium. In particular, the predicted shad- 
owing by the EPS parameterization is considerably larger 
than in the previous nuclear PDF's. The nuclear effects 
are different at central and at forward rapidities, so mea- 
suring the nuclear production of J/ip + 7 and T + 7 in 
these two regions could give additional insights about the 
correct nuclear gluon distribution. These complemen- 
tary observables to the inclusive production are useful 
to shed light on quarkonium production where the lack 
of a good benchmark, both for proton-proton collisions, 
where the actual mechanism of quarkonium production 
is not known, and proton-nucleus collisions where the rel- 
ative magnitude of the different effects (shadowing, nu- 
clear absorption) cannot be fixed by present experimental 
data. 
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